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' ABSTRACT 



> 



We measure the 2-point correlation function, ^ag> between galaxies and quasar absorption 
line systems at z < 1, using the dataset of Morris & Jannuzi (2006) on 16 lines-of-sight (LOS) 
with UV spectroscopy and galaxy multi-object spectroscopy (paper I). The measurements are 
made in 2D redshift space out to tt = 2Qh^^ Mpc (comoving) along the LOS and out to 
I a = 2h~^ Mpc projected; as a function of HI column density in the range iVni — 10^"^ ~ 

QQ , 10^^ cm~^ , also for CIV absorption systems, and as a function of galaxy spectral type. This 

extends the absorber-galaxy pair analysis of paper I. We find that the amplitude of the peak 
in ^AG 3.1 the smallest separations increases slowly as the lower limit on iVni is increased 
\0 '. from IQi^ to 10^^ cm - , and then jumps sharply (albeit with substantial uncertainties) for 

' TVhi > lO^^cni^^ . For CIV absorbers, the peak strength of ^ag is roughly comparable 

fH I to that of HI absorbers with iVni > 10^^'^ cm^^ , consistent with the finding that the CIV 

Q^' absorbers are associated with strong HI absorbers. 

We do not reproduce the differences reported by Chen et al. between ID ^ag measure- 
ments using galaxy sub-samples of different spectral types. However, the full impact on the 
measurements of systematic differences in our samples is hard to quantify. We compare the 
^ I observations with smoothed particle hydrodynamical (SPH) simulations and discover that in 

the observations ^ag is more concentrated to the smallest separations than in the simulations. 
The latter also display a 'finger of god' elongation of ^ag along the line of sight in red- 
shift space, which is absent from our data, but similar to that found by Ryan- Weber for the 
■ cross-correlation of quasar absorbers and Hl-emission-selected galaxies. The physical origin 

' of these 'fingers of god' is unclear and we thus highlight several possible areas for further 

investigation. 

Key words: galaxies - intergalactic medium, galaxies - quasars: absorption lines, galaxies - 
galaxies:haloes 
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1 INTRODUCTION 

We are investigating the nature and extent of the relationship be- 
tween gas, as probed by neutral hydrogen quasar absorption line 
systems, and the luminous baryonic matter traced by galaxies at 
cosmic epochs 2: < 1. In this paper we extend the work presented 
by Morris & Jannuzi (2006) (paper I), where we set out the moti- 
vation for the project, defined the 16 lines-of-sight (LOS) observa- 
tional dataset, and established the existence of an absorber-galaxy 
correlation using pair count analysis. We extend the analysis of this 
dataset to compute the absorber-galaxy 2-point correlation func- 
tion, ^AG, as a function of HI column density and compare the re- 
sults with other measurements and simulations. 

As discussed and referenced extensively in paper I, the history 
of work in this area has been characterised by a debate between two 
extreme positions: (i) that low-redshift Lya absorbers arise in the 



haloes of luminous galaxies; or (ii) that they are part of a filamen- 
tary network pervading the intergalactic medium (IGM), and are re- 
lated to galaxies only to the extent that they both trace the underly- 
ing dark matter distribution. This polarised picture is undoubtedly a 
simplification of a reality which now seems likely to require a mix- 
ture of the two components in proportions which vary as a function 
of galaxy type, absorber column density and cosmic epoch. A fur- 
ther complication arises from the expectation that at high redshift 
galaxies are expected to blow most of their star-forming material 
out into the IGM as "superwinds". Such outflows are hypothesised 
to pollute the IGM with metals (Aguirre et al. 2001), to prevent the 
formation of over-luminous galaxies at z = (Benson et al. 2003) 
and to account for the low-fraction of z = baryons locked in 
stars. The measurement of ^ag provides a method to characterise 
these feedback processes and to facilitate a direct comparison with 
simulations. At 2 ^ 2 — 3, attempts have been made to mea- 
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sure the impact of superwinds on the IGM by cross-correlating 
HI and CIV absorption systems with UV selected galaxies, but 
the early results were ambiguous (Adelberger et al. 2003). Sub- 
sequent measurements have clarified the situation (Adelberger et 
al. 2005), and reveal that strong Lya absorption is produced in 
intergalactic gas within lh~^ comoving Mpc of ~ 2/3 of such 
galaxies, and that strong CIV absorption (A^'civ 3> 10^^ cm~^ ) 
extends out to impact parameters of 40 kpc (proper) with a veloc- 
ity range Au > 260 kms^^ . Whether such superwind events have 
imprinted a relic signature on the galaxy-IGM relationship which 
persists to z < 1 is not known, but the measurements presented in 
the present paper may provide some contraints on this issue (see 
also Ranch et al 2005). 

In paper I we showed that absorbers and galaxies are corre- 
lated out to impact parameters of at least 1 .5 Mpc (physical) but that 
this clustering is weaker than the galaxy-galaxy clustering in the 
sample. The velocity differences between absorber-galaxy pairs 
with detected CIV absorption are typically smaller than those with 
HI absorption only, but variations with HI column density were not 
assessed. Several other works have cross-correlated absorption line 
catalogues with various low-z galaxy samples and reached broadly 
similar, but not identical, conclusions. Ryan-Weber (2006) cross- 
correlated galaxies selected in HI 21 cm emission with low column 
density (A^hi < 10^^ cm~^ ) Lya absorbers and found the clus- 
tering to be similar in strength to the clustering of the galaxies in 
the sample in redshift space, but somewhat stronger in real space 
on scales of l-\Qh^^ Mpc. At somewhat higher column density, 
Bouche et al. (2004) cross-correlated Mgll absorbers with lumi- 
nous red galaxies (LRGs) at 0.4 < z < 0.8, and over comov- 
ing scales of 0.05-13/i~^ Mpc found that the Mgll-LRG cross- 
correlation amplitude is ~ 30 per cent weaker than the LRG- 
LRG clustering. In a series of papers Prochaska and collaborators 
have investigated the clustering of z < 0.5 galaxies around ab- 
sorbers along the sightline to the quasar PKS 0405-123 (Prochaska 
et al. 2004; Chen et al. 2005; Prochaska et al. 2006). To the ex- 
tent that it is possible using a single sightline, they showed that 
emission-line galaxies and absorbers with A'hi > lO^** cm~^ clus- 
ter with a strength comparable to the galaxy-galaxy clustering, 
whilst the absorption-line galaxies in their sample are not cor- 
related with absorbers 0. They also showed that absorbers with 
A^Hi < 10^'^'® cm^^ are distributed more randomly with respect 
to galaxies, and that the strength of ^ag is insensitive to A^hi for 
^Qi3.6 ^ ^ -^Qi6.5 ^^-2 ^^^^ jjjig sightline, Williger 

et al. (2006) found that the absorber-galaxy correlation function is 
significant out to Au < 250kms~^ and grows with minimum 
absorber column density. Some of these conclusions can be tested 
with our multiple LOS dataset and that is our aim here. 

For completeness, we note that Chen et al. (2001) presented a 
sample of quasar lines of sight with HST FOS observations together 
with galaxy redshifts. As mentioned in Morris & Jannuzi (2006), 
we have chosen not to fold their data into our current analysis due 
to significant differences in survey strategy. 

Throughout the paper we adopt the cosmological parameters 
Q.M = 0.3, Q.A = 0.7, audi?,) = 100/ikms"^ Mpc"^ 



2 THE ABSORBER-GALAXY 2-POINT CORRELATION 
FUNCTION 

The absorber galaxy cross correlation function, .^ag, is defined 
through the expression for the conditional probability dP of find- 
ing a galaxy in a volume dV at position r2, given that there is an 
absorber at position ri : 



dP = ndV[l + aG(r2 -ri)], 



(1) 



where n is the unconditional galaxy density. 

The samples of absorption lines and galaxies which we will 
use to compute f ag are identical to those defined in paper I, where 
full details can be found. In brief, the dataset comprises 16 quasar 
fields (from the Hubble Space Telescope Quasar Absorption Line 
Key Project dataset; Jannuzi et al. 1998), whose LOS contain a total 
of 381 Lya absorption lines and 30 CIV systems (25 at 2; < 1), and 
a total of 685 galaxies with redshifts z < 1. As discussed in paper 
I, 49 of these galaxies have recession velocities < 500 km and 
are most likely stars, but their exact identification is irrelevant for 
this study because 98 per cent of the HI absorbers are in any case at 
z > 0.3. The number of galaxies in the redshift range 0.3 < 2 < 1 
which effectively overlap with the absorption line sample is 379. 

The following estimator, introduced by Davis & Pee- 
bles (1983), is used to compute ^ag as a function of projected sep- 
aration, a, and line-of-sight separation, n: 



iAG{o,Ti)= AG/ AR^l, 



(2) 



where AG is the number of real absorber-galaxy pairs in the 
region (ct — 5ct/2 : cr + 5a /2, n — Sn/2 : n + 5n/2) and AR is 
the suitably normalized number of absorber-random galaxy pairs at 
the same location. Following Bouche et al. (2004), the numbers of 
pairs in each of the quasar fields are normalized separately such that 
AR — '£AR''Ng/N^, where the summation is over the 16 fields 
and A'^^ and Nj. are the numbers of real and random galaxies in 
field i, respectively. All absorber-galaxy separations are computed 
in comoving redshift space and the absorbers are assumed to be 
located at the sky position of the quasar LOS. For the assumed flat 
cosmology, the comoving distance to an object at redshift z is given 
by: 



Xcm{z) = 



Ho Jo 



dz 



(3) 



and hence for a close absorber-galaxy pair separated by A2: in red- 
shift and by an angle A9 on the sky, a and n are given by: 



a = XcmA^ 



and 

dXcM . 

TV = ; AZ. 

dz 



(4) 



(5) 



Several crucial details of the computational method, namely 
the generation of the random galaxy catalogue, the binning and 
smoothing of ^ag, and the calculation of uncertainties, are dealt 
with in the next subsection. 



^ For conformity with earlier works, we retain the nomenclature of 
'emission-line' and 'absorption-line' galaxies to refer to galaxies whose 
spectra are dominated by emission line and absorption line features, re- 
spectively. The latter should not be confused with quasar absorption lines 
against which the galaxy populations are cross-correlated. 



2.1 Computational method 

Random galaxy catalogues 

Our method requires a separate random galaxy catalogue for each 
quasar field in order to reflect any field-to-field differences in the 
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galaxy selection function. The latter mainly arise due to the place- 
ment of slits and not from differences in spectroscopic depth. Fol- 
lowing Adelberger et al. (2003), we form such catalogues by re- 
placing each real galaxy with 100 random galaxies, whose redshifts 
are drawn at random from the galaxy redshift distribution of the 
other fields, i.e. the real and random galaxies have the same posi- 
tions on the sky, but different redshifts. This reflects the fact that 
galaxies in the sample must be located at the position of the slits. 



Spatial binning and smoothing o/^ag 

To reduce the effects of shot noise due to the small numbers of 
absorber-galaxy pairs on the measurement of ^ag, some consider- 
ation must be given to the form of the numerical coordinate grid, 
{a, tt), on which it is evaluated. There are two approaches. The first 
approach adopts a grid whose cell-size is large enough to include 
a sufficient number of absorber-galaxy pairs. An adaptively-sized 
mesh would be the best way to handle strong variations in the den- 
sity of absorber-galaxy pairs, but this would have the tendency to 
dilute the signal in regions where the paucity of pairs arises from a 
genuinely strong anticorrelation signal with ^ag — —1. 

An alternative approach is to bin the measurements much 
more finely (if at all) and then to smooth the resulting field to re- 
move the noise. Again, an adaptively-sized smoothing kernel could 
be the best option. The two underlying density fields, AG{a, n) 
and AR{a, n), must be smoothed first and then ratioed to obtain 
the smoothed ^ag, as in eqn. (2). 

In this work we present results obtained using both the 
above approaches: (i) on a uniform grid sampled as cr = 
0, 0.4, 0.8. . .2/1"^ Mpc, n = 2,4:, 6, ...20h'^ Mpc- referred to as 
the 'binned' results (the 2h~^ Mpc bin size in the vr direction re- 
flects the ~ 200 km s^^ spectral resolution; in the a direction, 
the field sizes dictate that galaxies are present out to maximum 
comoving impact parameters of ~ 2h^^ Mpc, which we have ar- 
bitrarily split into 5 equally spaced bins); (ii) on a finely-space 
grid with the underlying density fields smoothed by a gaussian of 
FWHM — 0.8 X Ah~^ Mpc, referred to as the 'smoothed' results. 
We did not experiment with any adaptively-sized grids or smooth- 
ing kernels. 



Uncertainties on ^ag 

Using the same method as Ryan- Weber (2006), we attempted to 
measure the uncertainties on ^ag using jackknife resampling in 
which the correlation function is computed = 16 times with one 
sightline removed each time. 

(6) 

A disadvantage of this method is that at locations where there 
are no absorber-galaxy pairs in the full sample (i.e. ^ag = — 1), 
all the individual = — 1 and hence cr^ = 0. This highlights an- 
other problem of assessing the significance of the measured ^ag 
values, namely that there is no obvious figure of merit analogous 
to a signal-to-noise ratio, i.e. ^ag is not a straightforward in- 
dication of a correlation's significance. An simpler alternative to 
characterise the significance of the observed correlation is to use 
Possion statistics, i.e. to calculate the probability of observing the 
real number of absorber-galaxy pairs (AG) in a given cell, given the 
expected number (AR). We will utilise both these approaches. 



3 RESULTS AND COMPARISON WITH OTHER WORKS 

In so far as the sizes of our absorber and galaxy samples permit, 
our goal is to measure ^ag as a function of both HI column density 
and galaxy spectral type, and to compare with other results in the 
literature. 



3.1 Dependence of absorber-galaxy cross-correlation 
function (^ag) on absorber HI column density 

The absorption line sample spans a range of HI column density 
from iVni = lO" - 10^^ cm"^ , with the bulk of the lines being 
at the lower end of this range, as shown by the distribution in Ta- 
ble 1. As discussed in paper I, the modest resolution of the HST 
Faint Object Spectrograph (FOS) (R=1300; 230 km s"^) means 
that the Lya absorption features were typically unresolved, so col- 
umn densities were inferred from the equivalent width assuming a 
Doppler 6-parameter of 30kms~^ . As shown in Fig. 1 of Wey- 
mann et al. (1998), the 4.5(7 lower detection limit on the rest 
equivalent width (REW) of absorption features in these quasars de- 
creases from O.2A al z = 0.3 to 0.1 A at 2; = 1. Considering that 
the overlap between the galaxy and absorber sample is skewed to- 
ward the lower end of this redshift range, the lower detection limit 
for the absorber sample is effectively REW ~ 0.2A, equating to 
Nui — 10^"^'^ cm~^ . We note, however, that the Key Project ab- 
sorption line sample is not strictly complete for any REW value, 
but has a detection limit which varies with both LOS and redshift 
(as discussed at length in paper I). 

In the first calculation of ^ag we use the entire absorption line 
sample, i.e. TVhi > 10^^ cm~^ . The results are shown in Fig. [T] 
where we also display the jackknife errors and the distributions 
of absorber-galaxy real (AG) and random (AR) pairs. The verti- 
cal banding in the distribution of random pairs results from an in- 
creased volume per cell at larger projected distances and a fall-off 
in the galaxy density at the largest radii in our sample. Observe that 
whilst the correlation function itself clearly peaks in the cell at the 
smallest separations, the peak in the absolute number of absorber- 
galaxy pairs is spatially offset, due to volume effects. This high- 
lights the value of basing inferences about statistical associations 
on ^AG rather than on the raw pair counts. The peak 'signal-to- 
noise' ratio, ^ag/o"{, is approximately 3. 

The evaluation of ^ag for various absorber sub-samples is 
presented in Fig. |2l with the same binning scheme as that used 
in Fig. [T] As the lower threshold on the HI column density is in- 
creased from 10" - lO^'^ cm , the strength of the central peak 
in ^AG also increases smoothly from ~ 4.5 to ~ 6.7. At A^hi > 
10^^ cm~^ the peak strength jumps dramatically to ~ 15 although 
the formal jackknife error is substantial. For absorber samples with 
a lower maximum HI column density, the peak in ^ag is of lower 
amplitude, ~ 3. For the absorbers with iVHi = 10^^ - lO^'^ cm"^ , 
the peak is spatially offset from the central position. These results 
are summarised in Table 2 where we list the peak values of ^ag, 
their jackknife uncertainties, and statistical significance assuming 
Poisson statistics. In Fig. [3] we show cross-correlation functions 
for the same HI column density intervals evaluated on a fine grid 
but with the underlying pair count fields heavily smoothed, as de- 
scribed in section 2.1. 

In Fig. [3] we also show the auto-correlation function of the 
galaxies in the sample. We note that whilst the peak value is com- 
parable to other measurements of the galaxy-galaxy 2-point cor- 
relation function on scales of l~2h^^ Mpc (e.g. from the 2DF; 
Hawkins et al. 2003), it is formally lower than the peak in ^ag for 
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Table 1. Distribution of HI absorber column density 



log A^Hi t Number of absorbers 



> 13 


379 






> 14 


241 






> 15 


161 






> 16 


91 






> 17 


39 






> 18 


16 






f Inferred assuming Doppler 


parameter 6 = 30 km s ^ . 




Table 2. Peak strength of ^ag shown in plots of Fig.|2] 




log Afni 


Peak 5ag 


Peak pair 


Peak 






real (random) counts* 


significance ** 


13-19 


4.4 ± 1.7 


13 (2.4) 


2 X 10-"^ 


14-19 


5.2 ± 2.5 


13 (2.1) 


4 X 10-'^ 


15-19 


6.0 ±3.8 


13 (1.9) 


1 X 10"'^ 


16-19 


6.7 ±5.2 


9 (1.2) 


4 X 10-*^ 


17-19 


15 ± 12 


7 (0.44) 


4 X 10"^ 


13-15 t 


3.0 ±2.0 


5(1.3) 


9 X 10-3 


15-17 


3.3± 1.1 


6(1.4) 


3 X 10-3 


CIV systems 


11 ± 9.9 


11 (0.94) 


1 X 10-* 



* Number of real (random) absorber-galaxy pairs in cell at ^ag peak 
** Poisson probability of observing > peak pair counts. These are included 
for completeness but substantially overestimate the significance level of 
any signal; the jackknife errors offer a more reliable significance measure, 
f Peak not in central bin 



the A'^Hi ~ lO^'^ — 10^^ cm^^ sample, which appears surprising. 
However, in view of the substantial uncertainty on the peak value of 
Cag for the latter subsample (see Table 2), any such differences in 
peak amplitude are not significant. We also note that other authors 
have compared the absorber-galaxy cross-correlation function to 
the galaxy-galaxy auto-correlation function to estimate the mass of 
the dark matter haloes in which the absorbers are embedded. Whilst 
this procedure is acceptable for Damped Lyman Alpha absorbers 
(which are likely to be galaxies) (e.g. Bouche et al. 2004), it is not 
appropriate for diffuse Lya absorbers below A''hi = 10^^ cm^^ . 
The latter are likely to arise in unbound cosmological filaments 
(e.g. Dave et al. 1999) so we do not pursue such calculations here. 
Indeed, Ryan- Weber (2006) used this method and found that weak 
Lya absorbers should arise in 10^* M© groups; this seems unphys- 
ical since gas in such groups would be quite hot, while the individ- 
ual galaxies' interstellar media do not have nearly the cross-section 
to reproduce the counting statistics of weaker Lya forest absorbers. 

For CIV absorbers, Fig.|2]shows that the peak strength of ^ag 
is comparable to that of HI absorbers with A^hi > 10^®'^ cm~^ . 
This is not surprising as in our dataset the median HI column den- 
sity of CIV absorbers is 10^^'^ cm^^ , while the median of the full 
Lyof absorber sample is 10^^'^ cm-^ . The median HI column den- 
sity of Lya absorbers with iVni > 10^*^'^ cm~'^ is 10^^ ''^ cm^^ . 

3.2 Dependence of ^ag on galaxy spectral type 

Chen et al. (2005) reported that the strength of the absorber-galaxy 
cross-correlation may be dependent on galaxy spectral type, in the 
sense that a strong correlation signal is found for emission-line 
galaxies but not for absorption-line galaxies. Here we investigate 
this issue by splitting our galaxy sample into two sub-samples, ac- 



cording to whether the redshift assigned in paper I was obtained by 
identification of a strong emission line or where no such emission 
line could be seen. The resulting emission-line and absorption-line 
galaxy samples contain 406 and 225 galaxies, respectively. Of the 
original full galaxy sample of 685, a total of 54 galaxies could not 
be classified and were not included in this analysis: 8 of them have 
spectra with an ambiguous classification, and the redshifts of 46 
galaxies were obtained from the literature without associated spec- 
tral information (see paper I for full details). The 225 absorption- 
line galaxies include most of the 49 very low redshift objects which 
were classified as probable stars in paper I, but this is not critical 
because such objects are too low in redshift to overlap with the ab- 
sorption line sample. 

These two galaxy samples were cross-correlated with the full 
HI absorption line catalogue (i.e. A'^ei > 10^"^ cm-'^ ) and the re- 
sults are displayed in Fig.|4] For the emission-line sample, the form 
of Cag is qualitatively similar to the result obtained with the full 
galaxy sample (Fig. [TJ, but the peak strength is now marginally 
(but not significantly) higher at ~ 5.8, compared with 4.4 ±1.7 for 
the full galaxy sample. With just the absorption-line galaxies, the 
central peak in ^ag is observed to be displaced from the central bin 
to 2-^ h^^ Mpc along the LOS, although we do not claim that this 
is necessarily a statistically significant finding, in view of the much 
smaller size of the absorption-line galaxy sample. A more direct 
comparison with the 1-D results of Chen et al. is presented in the 
next subsection. 



3.3 Comparison with other works 

Of the ^AG measurements mentioned in the Introduction, that by 
Chen et al. (2005) is the most similar to our work in its method- 
ology, as it is based on an optical magnitude-limited galaxy red- 
shift survey along the LOS to a single quasar (PKS 0405-123 at 
z — 0.5726). Here we make a comparison between our results 
and their 1-D measurements of ^ag as function of galaxy type. We 
begin with a summary of the Chen et al. dataset. Their galaxy sam- 
ple comprises 482 galaxies with i? < 20 mag and spectroscopic 
redshifts over a 1600 arcmin^ area within a 5h^^ Mpc (comov- 
ing) impact parameter of the quasar LOS; the redshift distribution 
has a broad peak at 2 ~ 0.2. For their ^ag analysis they restrict 
their attention to galaxies with impact parameter < lh~^ Mpc of 
the LOS, of which there are 61 (46 emission-line galaxies, and 
15 absorption-line galaxies; no details of the spectral classification 
method are given). The redshift distribution of this small impact 
parameter subsample is skewed towards lower redshifts compared 
with the entire sample, with a median redshift 2; ~ 0.1. They claim 
a galaxy redshift accuracy of Az — ±0.0002, or ~ 40kms^^ . 

The Chen et al. absorption line sample comprises 112 Lya 
systems at redshifts 0.01 < z < 0.557 with HI column den- 
sities logNm = 12.5 — 16.5 (80 per cent of the lines are at 
logNui < 13.6, the regime in which our absorption line sample is 
most incomplete). The bulk of the lines were derived from observa- 
tions with the HST Space Telescope Imaging Spectrograph (STIS) 
with spectral resolution of 6.7kms-^ . The latter is substantially 
higher than the ~ 230 km 3-^ resolution of the FOS spectra from 
which our absorption line sample was derived. Compared with our 
dataset, the Chen et al. sample can thus in principle probe ^ag on 
(i) smaller velocity scales; (ii) for weaker absorption lines ; (iii) 
with fewer complications arising from absorption line blending. 
The first two of these differences can be handled by comparing on 
appropriate Av scales and for sufficiently high TVhi, but the issue 
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Figure 1. For the complete sample of absorbers with A'^hi > 10^'' cm~^ , the figure shows the absorber-galaxy 2-point correlation function 
(top left), ^AG, determined using the real and random absorber-galaxy pairs counts (AG and AR, respectively, in eqn (2)) (lower left and 
right, respectively). Also shown is the uncertainty on ^ag derived using jackknife resampling, eqn. (6) (top right). 



of line blending in the FOS spectra is less straightforward and will 
be discussed later. 

Mindful of the above differences between the datasets, we re- 
strict our attention to a comparison with the Chen et al. ^ag mea- 
surements for absorption lines with logNm = 14.0 — 16.5. Chen 
et al. computed ^ag for galaxies within l/i~^ Mpc comoving of 
the LOS for the entire galaxy sample and its emission-line and 
absorption-line galaxy subsamples, computed out to LOS velocity 
separations of Av = 5000 kms"^ . Converting Av to comoving 
distance, we plot the Chen et al. measurements along with our re- 
sults out to LOS separations of 20h~^ Mpc in Fig.jS] Error bars on 
our data points were derived using jackknife resampling; Chen et 
al. give no details of their error bar computation. 

The plots in Fig. [5] demonstrate that we do not reproduce the 
high peak in ^ag at the smallest separations (comoving LOS sep- 
arations < 2h~^ Mpc) found by Chen et al. for the entire galaxy 
sample: our peak values are more than 2a apart. Possibly related to 
this, we also do not substantiate their finding of much stronger ^ag 
for the emission-line galaxies as compared with the absorption- 
line galaxies. Whilst our dataset is substantially larger than that 
of Chen et al. in terms of the numbers of galaxies (our study: 179 
galaxies with impact parameters < Ih^^ Mpc overlapping with the 



absorption-line redshift distribution; Chen et al: 61 galaxies) and 
absorbers (our study: 180; Chen et al.: 15), a much larger dataset 
will be needed to address more conclusively any dependence of 
^AG on galaxy type. It should also be borne in mind that the Chen 
et al. study is effectively at redshift 2 ~ 0.1 (in terms of maximal 
overlap between absorber and galaxy samples), whereas our study 
is at 2 ~ 0.5. Concerning the possible impact of different degrees 
of absorption-line blending in quasar spectra from the FOS and 
STIS, we refer to Penton, Stocke & Shull (2004): they smoothed 
their STIS/Goddard High Resolution Spectrograph (GHRS) spec- 
tra to FOS resolution and found a 15-25 per cent increase in the line 
density at REW> 0.24A (A^hi ^ 10^* cm"^ for b = 30 kms"^ ) 
due to the blending of weaker lines. Since Chen et al. found that 
these weaker absorption lines correlate much less strongly with 
galaxies, this alone would suppress our peak ^ag values by 25- 
35 per cent relative to their measurements, going some way to ex- 
plaining the differences in Fig. [5] However, the full impact on the 
^AG measurements of differences between our respective galaxy 
and absorption-line samples merits further investigation. Finally we 
note that the 1-D ^ag in Fig.|5]all exhibit possible minima at LOS 
separations ~ 10^~^ Mpc, similar to the feature in the smoothed 
2-D plots of Fig.|3] the physical origin of this feature is not known. 
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Figure 2. ^ag evaluated in 2-D redshift space for various ranges in HI column density (as indicated) and for CIV systems. Note the variations 
in the colour bar scales, and refer to the estimated uncertainties in Table 2 when comparing the peak ^ag values of the subpanels. 
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Figure 3. ^ag evaluated on a finer grid tlian sliown in Fig. [2| and smoothed as described in section 2.1. Results are presented for various 
intervals of HI absorber column density (as indicated) and the galaxy-galaxy clustering is also shown. Note the variations in greyscale value, 
and refer to the estimated uncertainties in Table 2 when comparing the peak ^ag values of the subpanels. 
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Figure 4. ^ag evaluated for the emission and absorption-line galaxy sub-samples using the full sample of absorption lines with iVni > 
10^^ cm~^ , as described in section 3.2 



Ryan- Weber (2006) cross-correlated quasar absorption line 
systems with 12.4 < logNui < 14.8 against an Hl-selected galaxy 
catalogue from the HIPASS survey. As in the Chen et al. absorption 
line sample, this absorption line sample is based on HST STIS and 
GHRS spectra and thus provides considerably higher velocity res- 
olution and weak line sensitivity compared with our sample. The 
galaxy sample is biased towards gas-rich galaxies, including gas- 
rich dwarf and low surface brightness galaxies often absent from 
magnitude-limited optical catalogues such as ours. These differ- 
ences, together with the lower redshift of her measurements (which 
are effectively at 2; ~ 0), might account for the extreme LOS elon- 
gation seen in her ^;ac; measurements - out to lOh^^ Mpc comov- 
ing - which we do not reproduce. Indeed, morphological analysis 
of the optical covmterparts to the 1000 brightest HIPASS sources 
confirms that the majority are late-tj^je gas-rich galaxies (Zwaan 
et al. 2003). The origins of this LOS elongation are explored fur- 
ther in section 4.2 below with the help of SPH simulations. In 
seeking to compare our ^ag results with those of Ryan- Weber, 
an important first step is to compare the ranges of galaxy mass 
probed by the two surveys. Ryan- Weber states that the geometric 
mean HI mass of galaxies contributing to absorber-galaxy pairs is 
log{Mm/ Mq) = 8.8h~^, with a corresponding dark halo mass 
of log(M/ Mq) = 11.0/i~^ (as inferred from modelling by Mo 
et al. 2005). In contrast, our galaxy sample was selected by R-band 
imaging and follow-up spectroscopy and is thus expected to contain 
a higher fraction of early-type galaxies. In paper I (section 2.3.5) 
we attempted to place bounds on the rest-frame B-band absolute 
magnitudes of our galaxy sample: confining our attention to galax- 
ies at redshifts 0.35 < z < 0.8 which overlap with the absorption 
line sample, we find that the vast majority of these galaxies have 
— 22 < Mb < —18 mag, with a median Mb — —20.5 mag. 
With reference to models of the conditional probability distribution 
of halo masses (e.g. Cooray 2005 and references therein), we ex- 
pect galaxies at this median luminosity to be hosted by dark matter 
haloes of mass - 10^^"^^/i"^ Mq. 

We note that not only are our galaxies hosted by more mas- 
sive haloes than those of Ryan- Weber, but that the two samples fall 
either side of a critical halo mass, at which the nature of gas ac- 
cretion onto galaxies is thought to change. Using SPH simulations, 
Keres et al. (2005) find that galaxies hosted by haloes with Mhaio < 
10^^-* M0 (or, equivalently, baryonic mass Mgai < 10^°-^ M©) 
acquire most of their gas via 'cold mode' accretion in which gas 



at T < lO^K is efficiently drawn in from large distances along 
filaments. For more massive galaxies, the dominant mechanism is 
quasi-spherical accretion of hot shocked gas at the virial temper- 
ature of the halo, more akin to the classical picture outlined by 
White & Rees (1978). Moreover, this transition in accretion mode 
is a strong function of halo mass over the range of interest: from 
A4aio = 10" M© to 10^^ M0, the fraction of gas accreted by 
the 'cold mode' accretion varies from > 0.9 to < 0.1, virtually 
independent of redshift. Similar results have been fovmd by Dekel 
& Bimboim (2006). We do not wish to claim that this switch from 
cold mode filamentary accretion in Ryan- Weber's HIPASS galaxies 
to hot mode accretion in our galaxies can wholly explain the differ- 
ent degrees of LOS elongation in ^ag measurements using the two 
samples, but it may be a useful signpost to fruitful areas of future 
investigation. 



4 COMPARISON WITH SPH SIMULATIONS 

In this section we use cosmological simulations to generate syn- 
thetic 2D absorber-galaxy correlation functions to compare with 
the observations, and to gain insight into the physical origin of fea- 
tures such as the 'finger of god' which is seen in the results of 
Ryan- Weber (2006). 

We ran a simulation using Gadget-2 (Springel & Hemquist 
2003) with improvements as described in Oppenheimer & Dave 
(2006), in a cubic comoving volume of 64/i~^Mpc on a side. We 
employed the momentum-driven outflow model 'vzw' as described 
in Oppenheimer & Dave , and a cosmology concordant with 3rd- 
year WMAP results (Spergel et al. 2006): ilm = 0.26, ^a = 0.74, 
Ho = 71 km/s/Mpc, as = 0.75, and Clb = 0.044. The force reso- 
lution is 5/i.~^kpc, and the mass per gas particle is 2.72 x 10* M©. 
The galaxy finding procedure is described in Finlator et al. (2006), 
using a minimum resolved galaxy stellar mass of 8.7 X 10'' Mq. 

The simulation volume was pierced by 1000 randomly-located 
sightlines to simulate quasar spectra, and absorption line profiles 
were fitted using Auto VP (Dave et al. 1997). Absorber-galaxy pairs 
were measured and tabulated from these spectra for LOS sepa- 
rations At) < 1600 km and projected separations of Ar < 
16 Mpc (physical). Following conversion of pair separations into 
comoving redshift space, these catalogues were used to create 2-D 
absorber-galaxy correlation functions. Two such absorber-galaxy 
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Figure 5. A comparison of our ^ag measurements (thick lines) 
with those of Chen et al. (2005) (dashed lines) for different 
galaxy subsamples and absorption lines with 10^* < A'^hi < 
lO^'''^ cm~^ (upper panel: all galaxies: middle panel: emission- 
line galaxies; lower panel: absorption-line galaxies). 1 + ^ag is 
plotted as a function of comoving line-of-sight distance for galax- 
ies within Ih^^ Mpc of the line of sight. Error bars on our data 
points have been displaced in the x-direction for clarity. 



pair catalogues were created in order to compare with our observa- 
tional results and those of Ryan- Weber, as we now describe. 



4.1 Timing the simulations to compare with the 
observational results 

Firstly, to compare with the observational results presented in sec- 
tion 3 of this paper, the simulation output was extracted at 2 = 0.5 
(giving a total redshift pathlength of Az = 26.4). The pixel scale 
of the quasar spectra was set to 70kms^^ and convolved with a 
Gaussian of 230 km s^^ FWHM, in order to match the resolution 
of the //5r Faint Object Spectrograph. HI absorption lines were se- 
lected with rest equivalent width REW > 0.056A, corresponding 
to column densities A'^hi > 10^'' cm^^ for b — 30kms^^ . The 
strength of the ionizing background radiation was increased by a 
factor of 4 over the Haardt & Madau (1996) value in order to nor- 
malise the dN/dz value for REW > 0.24A lines to the HST Ab- 
sorption Line Key Project results (Weymann et al. 1998). The en- 
tire resolved galaxy population was used for the comparison, with a 
galaxy density of 0.01 145 per cubic h^^ Mpc comoving above the 
8.7 X 10^ Mq baryonic mass limit. As shown by, for example, the 
Deep Evolutionary Exploratory Probe 2 (DEEP2) Galaxy Redshift 
Survey (Willmer et al. 2006), this is comparable to the space den- 
sity of galaxies at this redshift in an absolute range which matches 
our observed galaxy sample (—22 < Mb < — 18 mag; see section 
3). 

The second absorber-galaxy pair catalogue was extracted from 
the simulations at redshift z = (giving a total redshift path- 
length of Az = 21.3) to compare with the Ryan-Weber results. 
The spectral resolution of the artifical quasar sightlines was set 
to 8 kms^^ to mimic the higher resolution HST STIS and GHRS 
spectra, and HI absorption lines were extracted for REW=0.0126- 
O.45A, corresponding to 10^^ * < A^hi < lO" ** cm"^ for 
b = 30kms~^ . Mimicing the HIPASS selection of galaxies in 
HI emission by means of a simple cut in baryonic mass is a crude 
approximation, but it is the only tool at our disposal with the current 
simulation. Somewhat arbitrarily, we use an upper baryonic mass 
cut-off for galaxies of 2 x 10^" M0 . To put this figure in context, we 
refer to some results from Durham semi-analytic models of galaxy 
formation (Cedric Lacey, private communication; for the latest im- 
plementation of such models see Bower et al. 2006). Such models 
track the behaviour of gas and stars in dark matter haloes with sim- 
ple physical prescriptions. At redshift z — 0, dark halo masses of 
10", 10^^, 10" and lO^^h'^ M© contain galaxies of total bary- 
onic mass (hot gas, cold gas and stars) of 10^°, 8.6 x 10^°, 10^^ 
and 10^^ M0, respectively. The corresponding cold gas masses are 
1.7 X 10®, 6.7 X 10®, 3.3 x 10^" and 2.7 x lO" Mq, respectively. 
A baryonic mass cut off of 2 x 10^° M0 corresponds to a cold gas 
massof4.2x 10® M© and a dark halo mass of 2.5x 10"/i"^Mq. 
As stated in section 3, the HIPASS galaxies contributing to Ryan- 
Weber's absorber-galaxy pairs have mean HI and dark halo masses 
of 10** **/^^ and IO^/ZiMq, respectively. The space density of 
these simulated galaxies is 0.006107 per cubic /i~^Mpc comov- 



4.2 Results: the simulated ^ag and the origin of the 'finger of 
god' 

The two-dimensional ^ag which result from the simulations de- 
scribed in section 4.1 are displayed in Fig. |6] Plots with two 
types of binning/smoothing are presented, firstly binned coarsely 
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Figure 6. ^ag for the two simulations described in section 4.1, each plotted in two different ways, (a) The simulation with STIS-resolution 
sightlines, cross-correlating 10^^'* < A^hi < 10^* ** cm^^ absorbers with galaxies of baryonic mass < 2 x 10^" M0 at z = 0. (b) 
As in (a), but binned/smoothed as in Ryan-Weber (2006). (c) The simulation with FOS-resolution sightlines, cross-correlating TVhi > 
10^"^ cm~^ absorbers at z = 0.5 with the entire resolved galaxy population; (d) as in (c) but binned/smoothed as in Ryan-Weber (2006). 



to match our results in Fig. |2] (although the simulated plots ex- 
tend out to slightly smaller LOS distances), and secondly using the 
technique of Ryan- Weber. The latter evaluates ^ag on a fine grid 
(0.1 X Q.lh~^ Mpc) and then smoothes the resulting ^ag by 9 x 9 
pixels. Note that this differs from the smoothing procedure we used 
to generate Fig. [3] where the underlying pair counts were smoothed 
before computing ^ag • 

Comparing the simulations of Fig. [6] with our observational 
results in Fig. [2] and with the results of Ryan- Weber, we notice 
that both sets of observed ^ag are much more strongly peaked 
at the smallest separations than the simulations, i.e. we observe 
much more gas at small velocity offset and/or close to the galax- 
ies than present in the simulations. Redshift measurement errors 
in the observations would tend to weaken such peaks, so this can- 
not be the origin of the discrepancy. A related issue is the asym- 
metric elongation of the simulated ^ag in redshift space along the 
LOS, which is qualitatively similar to that found by Ryan- Weber 
but which appears to be absent from our measurements. Such a 
feature is analogous to that seen in the autocorrelation function of 
galaxies (see e.g. Hawkins et al. 2003), where it arises from de- 



partures from Hubble flow due to virialised motions in groups and 
clusters. Commonly referred to as the 'finger of god', the effect was 
seen in ^ag in an earlier generation of cosmological simulations 
(e.g. Dave et al. 1999) but its physical interpretation is not as clear 
as in the galaxy-galaxy case. Ryan- Weber (2006) interpreted it as 
arising from the draining of gas from low-density regions into col- 
lapsed structures. In a subsequent paper we will address the role of 
galaxy outflows in causing the 'finger of god', and also investigate 
the dependence of ^ag on galaxy mass and absorption line column 
density. A more ambitious goal is to simulate the full observational 
dataset, including the selection of the galaxies in a somewhat more 
realistic manner than is possible with the existing cuts in total bary- 
onic mass. 

In closing this section, we comment on another possible ex- 
planation for the 'finger of god' which merits further investigation. 
It may be a geometrical effect which leads to a bias towards observ- 
ing large numbers of galaxy-absorber pairs at small projected sky 
separations when a LOS is oriented directly down an IGM filament 
with clumpy HI absorption. At redshifts z < 2, simulations show 
that absorbing gas with TVhi < 10^''' cm~^ resides in a filamentary 
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network (e.g. Dave et al. 1999), whereas at higher redshifts such 
gas is distributed more diffusely. This basic structure evolves little 
over Q < z < 2 but absorbing material of a given column density 
becomes increasingly clumpy towards lower redshifts. This might 
partially account for the increased LOS elongation measured by 
Ryan-Weber at z ~ compared with our measurements at 2: ~ 0.5, 
but further investigations are needed. If only a very small propor- 
tion of sightlines are fortuitously aligned with filaments, it is much 
more likely that such chance alignments occur in the simulations 
(1000 LOS; A2 ~ 25) which have a much higher redshift path- 
length than our dataset (16 LOS; ~ 3) or that of Ryan-Weber 
(27 LOS; ~ 1). 



5 CONCLUSIONS 

In summary, our measurements show that the strength of the corre- 
lation function between HI quasar absorption systems and galaxies 
increases with HI column density, certainly for comoving separa- 
tions < 0.4/i~^Mpc (projected) and < 2h~^M^c (line-of-sight) 
where the signahnoise of our measurements is greatest. The peak 
strength of the absorber-galaxy clustering increases only slowly 
as the lower TVhi limit is increased from 10^^ to 10^^ cm~^ , but 
then jumps considerably (albeit with substantial uncertainty) above 
10^^ cm~^ . If real, such a jump may reflect the transition from ab- 
sorption within the filamentary IGM to absorption within the haloes 
of individual galaxies. The peak strength of the CIV absorber- 
galaxy cross-correlation is comparable to that of the HI absorber- 
galaxy signal for iVni > lO^*^'^ cm"^ . 

Our results are broadly consistent with the finding of Chen et 
al. (2005), based on the sightline to the quasar PKS 0405-123, that 
the absorber-galaxy cross-correlation strength is insensitive to HI 
column density in the range lO" '' < iVni < lO^'''^ cm"^ , al- 
though Table 2 shows evidence for a weak increase with A^hi . We 
find tentative but inconclusive evidence for differences in the 2-D 
absorber-galaxy cross correlation function for our subsamples of 
emission-line and absorption-line galaxies, but we do not substan- 
tiate the claim of Chen et al. that the 1-D cross-correlation signal 
is much higher for emission-line galaxies than for absorption-line 
galaxies. This discrepancy between our results and those of Chen et 
al. is not fully understood, but it may be partially ascribed to the ef- 
fects of line blending in the FOS spectra compared with the higher 
resolution STIS quasar spectra. 

Comparison was also made with the results of Ryan- 
Weber (2006) on the cross-correlation of Hl-selected galaxies with 
quasar absorbers with 10^^ * < iVni < lO" * cm , who find 
a marked 'finger of god' elongation out to IQh^^ Mpc along the 
sight. Such a feature appears to be absent from our measure- 
ments, but we demonstrated that there are significant differences 
in the characteristic mass of our respective galaxy samples which 
may account for such differences. The dark matter halo masses 
of the two samples fall either side of a critical dark halo mass of 
lO^i-* M0 (with the Ryan-Weber sample being lower in mass), at 
which simulations show that the nature of the dominant mode of 
gas accretion onto galaxies changes fundamentally from cold fila- 
mentary accretion in low mass haloes to hot quasi-spherical accre- 
tion (Keres et al. 2005). 

We also constructed synthetic absorber-galaxy 2-point corre- 
lation functions using output from a 64/i^^ Mpc box SPH simula- 
tion aimed at reproducing our own observational results and those 
of Ryan- Weber (2006). The simulations revealed somewhat less 
gas at the smallest separations from galaxies than present in the 



observations and a strong 'finger of god' elongation of ^ag along 
the LOS. These discrepancies may originate in true deficiencies in 
the physics of the SPH simulations, or in a failure to accurately 
mimic the observed galaxy selection functions, or through some 
more subtle selection effect. 

For future work, we suggest that effort be invested to under- 
stand the dominant contribution to the 'finger of god' elongation 
in ^AG, and the nature of the apparent transition in ^ag between 
iVni = 10^® and 10^^ cm . On the theoretical side, comparison 
should be made with a variety SPH simulations matching the obser- 
vational selection function more closely. On the observational side, 
it would also be valuable to conduct surveys for higher-redshift 
galaxies around existing HST Quasar Absorption Line Key Project 
sightlines to ensure a better match between the redshift distribu- 
tions of the absorption lines and the galaxies, and around higher- 
resolution HST STIS quasar sightlines. It would also be valuable to 
survey some of these sightlines with large-format IFUs to enhance 
sensitivity to emission-line galaxies which may be preferentially 
missed by our existing broad-band optical selection techniques. 
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